Density functional methods for polymers: a coil-globule transition case study 
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We consider a free energy functional on the monomer density function that is suitable for the study 
of coil-globule transition. We demonstrate, with explicitly stated assumptions, why the entropic 
contribution is in the form of the Kullback-Leibler distance, and that the energy contribution is 
given by two-body and three-body terms. We then solve for the free energy analytically on a set of 
trial density functions, and reproduce de Gennes' classical theory on polymer coil-globule transition. 
We then discuss how our formalism can be applied to study polymer dynamics from the perspective 
of dynamical density function theory. 

PACS numbers: 82.35.Lr, 71.15.Mb 



I. INTRODUCTION 

Polymer coil-globule transition at equilibrium has im- 
portant biological and technological importance and has 
thus received its well-deserved attention P, H, H, 0, ID, d, 
0, H, [H El- Here, we attempt to devise a simple for- 
malism that treats the free energy as a functional on the 
monomer density. The idea of having a free energy func- 
tional is of course not new 0, 0] , but we believe that our 
contribution distinguishes itself by its simplicity, and by 
its emphasis on the entropic ground state. Specifically, 
we consider a free energy functional of the form: 



T[p{r)] = / dr/(p(r)) 



(1) 



where / is some function, p(r) is the normalized 
monomer density function (i.e., J drp(r) = 1) and r 
is the Monomer-to-Center-of-Mass (MCM) distance (c.f. 
Fig. [T|) In constructing the free energy functional, we ar- 
gue that the entropic contribution should correspond to 
the Kullback-Leibler distance [TJ E3 between p(r) and 
the entropic ground state; and the energy contribution 
should consist of two-body term and three-body terms. 
We then introduce a set of trial density functions that 
are suitable for the coil-globule transition, and solve for 
the corresponding free energy analytical. We find that 
our approach is equivalent to de Gennes' classical the- 
ory on polymer coil-globule transition. We then discuss 
how our work connects to the excitin g de velo pm ent in 
dynamical density functional method jl3l . [lH HE) , and 
derive a second-order differential equation governing the 
dynamics of p(r) that incorporates the entropic effect, 
monomer-monomer interaction, brownian motion of the 
monomer-monomer interaction. 



II. THE FREE ENERGY FUNCTIONAL 

Our aim is to construct a free energy functional that 
depends on the MCM distance distribution, p{r). We 
firstly discuss the entropic part, and we start by writing 



FIG. 1: The distance r is the measured from the monomer to 
the center of mass of the polymer as depicted schematically 
below. 




the entropy term as: 

S = - ^pdogpi , 



(2) 



where the sum is over all available states. In a 3D regular 
lattice, a N segments phantom chain would have a total 
of Q = 6 N states. Let us denote the corresponding MCM 
density function for each of the path by 6i(r), 1 < i < CI. 
The optimal MCM density function is then: 



Po 



(3) 



We now assume that Qi are slowly varying in i and adopt 
the continuous notation: 



Po(r) = / dsO(s,r) 



(4) 



where s = i/Q and is in the range [0, 1]. For any density 
function p(r) other than po( r ) ; we have 



p(r) = Q / dsq(s)6(s,r) 



(5) 



where J dsq(s) = Q , and q(s) is not uniform as p ^ po 
by assumption. The corresponding entropy for p is 



fi / dsq(s)log(q(s)) 
Jo 



(6) 
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To construct an entropic functional S[p] of the form: 

S[p(r)] = [ drf(p(r)) (7) 



where / is some function, we need to go from the con- 
figuration enumeration index s to the MCM distance r 
(c.f. 0). As far as we know, there is no known analytical 
expression relating s and r. To make progress, we will 
now make a series of approximations. Firstly, we assume 
that there exists an inverse, 9, in Eq. [5]in the sense that: 



q(s) = / dr9(s,r)p(r) 



(8) 



We then assume that 9(s, r) is of the form 6(s)8(£(s) — r) 
where £ is a map that takes the configuration enumera- 
tion index s to the MCM distance r, and 8(.) is the Dirac 
delta function. The assumption is equivalent to saying 
that q(s) is determined by the value of p at the point 
r = £(s) alone. With this assumption, we have 



q (s) = e(s)p(m) 



which gives: 



S[ P ] = n / d s 0( s )p(£( S ))log[0( S )p(£( S ))] 
Jo 

Since we know from Eq. [4] that 

1 = 9(*)po(C(*)) - 
we can solve for 9(s) and obtain the following: 



Po(£(s)) 



p(£(*)) 



npo(f(*)) 



(9) 



(10) 



(11) 



(12) 



Let us now consider what the map £ should be. The 
interchange of the enumeration index and the spatial pa- 
rameter can be seen as a way to weight the sum in the 
above integral differently as the MCM distance varies. 
It is intuitive to set the weight according to the magni- 
tude of the optimal density function, po{r). Namely, the 
difference between p and po at r is given more empha- 
sis when po(r) is large. This intuition suggests that we 
set ds = pa(r)dr. In other words, the map £ is defined 
implicitly by the following equation: 



■s - / dr'p (r') 
'o 



(13) 



Now with this particular map, we have for the entropic 
functional the following expression: 



/>oo 

S[p] = - log n + / drp(r) log 
Jo 



p( r ) 
Po{r) 



(14) 



where we recall that log 51 = AHog(6). Besides the 
constant term, the above expression is identical to the 



Kullback-Leibler distance, or the relative entropy, be- 
tween p and pq. We note that the Kullback-Leibler dis- 
tance is a well studied measure in the context of informa- 
tion theory [ll|, Q3 and we have also recently emp loyed 
it in the context of sampling and optimization [161 ] . 

As a functional in p, a desirable property is that one 
can obtain the optimal density by simply considering the 
zeros of the corresponding functional derivative. Now, 



op \Po 



1 



Setting the above expression to zero suggests that: 



P = c Vo 



(15) 



(16) 



which is contrary to our expectation that p should be 
po exactly. This is a manifestation of the fact that the 
functional in Eq. [14] alone is not probability preserving. 
Indeed, given a system with M possible states, a better 
expression for the system's entropy is: 



5' 



M-l 

■E 

i=l 



Pi log(Pi) 



M-l 

E 

i=l 



Pi log 1 



M-l 

E; 

i=l 



(17) 

where the unity of the sum of a probability distribution 
is explicitly accounted for. With large M the second 
term in the above expression is small and as such the 

leading order term is ^1 — J^iLi 1 P^j- I n other words, 

the entropic functional can be expressed as: 



s[ P ] = -iogn + i 



drp(r) 



log 



P( r ) 
Po(r) 



(18) 

With this modification, it can be checked that po does 
correspond to the state where the functional derivative 
vanishes. Note that for a probability distribution p(r), 
the above modification is superfluous, but the manip- 
ulation may be necessary when considering dynamical 
effects. In summary, Eq. [T5] constitutes the entropic 
functional that we will use here. 

We now discuss the energy contribution in the 
free energy functional. We assume nearest neighbor 
interactions in our lattice model and we let rj be the size 
of the monomer. Consider a concentric spherical slice at 
position r of thickness dr. The number of monomers in 
this shell is: 

Np{r)dr , (19) 
and the number of sites available in this shell is: 
Airr 2 dr 



(20) 



Therefore, the probability of any particular site being 
occupied is: 



P(r) 



Np(r)if 
Aitr 2 



(21) 
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Adopting a mean-field perspective and ignoring the in- 
teraction between different shells for the time being, we 
have the following energy functional: 



EM = iefdr 



Attt 2 ( Np(r)rf 



rf 



N 2 p{r) 2 rj 



Airr 2 

2^,3 



(22) 
(23) 



where the factor 4 in front corresponds to the maximum 
number of neighbors allowed in our slice of volume. 

The interactions between the shell at position r and 
the shells at positions at positions r + rj and r — rj can be 
written as: 



e[p(r)p(r + rj) + p(r - T))p(r)] 
ep(r)[p{r + rj) + p(r — rj) — 2p(r) 

„2 // 



Mr) 



= ep{r)[r, 2 p"(r) + 2p(r)} 



(24) 
(25) 
(26) 



where p"{r) = d 2 p/dr 2 and rj is assumed to be small 
enough that the approximation by the differential oper- 
ator is possible. Writing it as an energy functional, this 
shell-shell interaction term becomes: 



E-2[p] = ~j drNp(r) 



dr 



( Np(r)r] 3 
\ 4irr 2 



N p{r)rj 
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2irr 2 



(27) 

where the factor 1/2 appears to take care of the double 
counting. We restrict ourselves to the case where rj is 
small, hence we will ignore the first term in the squared 
brackets in Eq. [27J which is of order rf . The overall 
energy contribution to the free energy functional is then 
simply: 



E[p] = | J dr 



N 2 p{r) 2 rf 



(28) 



We note that we could have arrived at the above func- 
tional by assuming that p(r) is slowly varying and so that 
the p"{r) term can be ignored; but this route would not 
allow us to see what the magnitude of error in our ap- 
proximation is. Our consideration also allows for the pos- 
sibility of investigating how the situation would change 
with large rj, in which case the energy term will be a 
functional dependent on both p and p" . 

As seen by inspecting Eq. [2H1 a two-body attractive 
term alone will drive the polymer to the complete collapse 
state at any finite temperature. This can be remedied by 
the volume exclusion effect. To account for such an effect, 
we assume that the collapse state corresponds to a ball 
with a uniform density, i.e., p(r) cx r 2 . This suggests 
that the energy functional should be of the form: 



E[p] = j dr 



Ap 2 



Bp 3 



(29) 



where A and B are arbitrary parameters. The second 
term above clearly corresponds to a three-body repul- 
sive term as usually introduced by hand in the literature 



FIG. 2: The circles denotes the shapes of po(r) correspond- 
ing to a phantom chain with unit length segments for various 
TV obtained from numerical simulations. The solid lines de- 
pict the approximation by </J 7 (jv)(r). For each TV, 7 is found 
by equating the mean of 4>~t to the mean of po- It is found 
that 7 = 0.0063,0.0031,0.0016 for N = 500,1000,2000 re- 
spectively. In other words, viewing 7 as a function of TV, we 
have 7(/V) x TV ~ n. 




(e.g., see [1JJ)- In our case, a three-body repulsive term 
arises naturally by fixing the ground state of the energy 
functional. 

In summary, by introducing two new parameters, w\ 
and u>2, our energy functional becomes: 



E[p] = Jdr 



N 2 p 2 



(wi - e) 



W2- 



N 3 p 3 



(30) 



where the constants in the attractive two-body term are 
absorbed into e. The optimizing density function for 
Eq. (3D] now corresponds to a uniform ball with density 
2(e - wi)/3u;2TV- 

Putting the entropic and energy terms together, 
the full free energy functional is: 



T[p] = -Tlogfi + T+ J drTp 

N 2 p 2 



dr 



(wi - ()- 



108 [fJ- 1 

N 3 p 3 
-W2—7- 



(31) 



III. VARIATIONAL METHOD 

We now apply the free energy functional to investigate 
polymer coil-globule transition. Looking at Eq. [5T1 we 
see that we firstly need to get a handle on po(r). We 
are unaware of any analytical solution for Pq(t). Here we 
adopt the simple assumption that for the a chain with TV 
segments of length l p , 



Po(r) 



4 7 3/2 
7TV2 \l 



(32) 
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where the corresponding radius of gyration is l p y^3/2~/ 
and as such, 



Rg 



(33) 



Fig. [2] illustrates that the approximation is in reasonable 
agreement with the density function obtained by numer- 
ical means. In particular, we find that 'y(N) x ./V ~ tt. 

The above consideration motivates us to introduce the 
set of trial density function of the form: 



Indeed, these maps could have been anticipated by di- 
mensional analysis on Eq. [33J and Eq. ??. These corre- 
spondences show that our variational method on the free 
energy functional are equivalent to de Gennes' theory 
on coil-globule transition, with except that the argument 
Rg in our approach is an index for the density function 
4>2/3R2,( r )- Namely, unlike de Gennes' approach, the free 
energy is no longer defined by the EED. Furthermore, 
our approach did not rely assuming that the radius of 
gyration is equal to the EED in the derivation. 



Mr) 



-gr 2 /il 



(34) 



IV. CONCLUSION AND OUTLOOKS 



where g is treated as the variational parameter. Since 
the radius of gyration corresponding to 4> g {r) is l p y/3/2g, 
which goes to zero as g goes to infinity, this set of density 
functions is also suitable to describe the collapse state. 
With that said, we note that these trial functions do not 
represent well the minimal energy configuration because 
a uniformly dense ball has sharp density cutoff at the 
boundary while our set of trial functions have stretched 
tails, but we believe that this discrepancy would not be 
very important when the polymer is away from the com- 
pletely collapse state. 

Besides being capable of describing both the free and 
collapse states, our trial functions also have the virtue 
of allowing us to compute the free energy in Eq. [3T] an- 
alytically. Solving for all the Gaussian integrals in the 
functional and writing the free energy as a function of g, 
we obtain: 



Ha) 



T 













+ (w 1 ~e)g^ 2 N 2 J- + 



3(g-7) 

7T 



w 2 - 



logfl + 1 

l6g 3 N 3 
3 3 / 2 tt ' 



(35) 



In this paper, we have introduced a free energy func- 
tional on the MCM density function, p(r), and applied it 
to the study of polymer coil-globule transition. In con- 
structing the free energy functional, we argued that the 
Kullback-Leibler distance measure is, under clearly delin- 
eated assumptions, a suitable functional for the entropic 
part. We then solve for the free energy exactly on a set 
of trial density functions and found that our approach 
is equivalent to de Gennes' classical theory on polymer 
coil-globule transition. 

In terms of outlook, we believe that our work makes a 
connection to the important development in dynamical 
density functional method [H, 0, [l5[ by allowing for the 
possibility of studying polymers dynamics with our for- 
malism. For instance, for any polymer out of equilibrium 
at t = 0, according to the dynamical density functional 
theory, the temporal evolution of the MCM density func- 
tion may be described by 



<3pM) = _± 

^ dt 



p(r,i)V 



5 P (r,t) 



(40) 



Recall that the radius of gyration, Rg, corresponding to 
4> g is \j2>j1g, we rewrite the free energy in terms of Rg 
and keep only terms that depend on Rq: 

jTR% , ^N 2 



HRg) 



N 3 

3T log(ifc) + + (^ - e ) _ + W2 _ 

G 1 



G 

(36) 

As a comparison, we reproduce the free energy formula 
as deduced originally by de Gennes below (c.f. Eq. 2.42 
in (3): 



HRc 



-3Tlog(i? G )- 



9TR 2 G 3{Tv - e) N 2 v 2 T N 3 
"lzfi\T + 4 R^ + 4 R% 

(37) 

where v is the volume of the monomer and e is the attrac- 
tive energy. By inspection, we see that the two formula 
are equivalent with the following mapping: 



7 



9 

w 

3e 
-> — 
4 



w 2 



3vT 
v 2 T 



(38) 
(39) 



where C is the friction constant. For our free energy func- 
tional, it leads to the following second-order differential 
equation: 



^dp = d_ 
dt dr 



4 f Ti °g- 

or \ po 



2(wi - e)N 2 p 3w 2 N 3 p 2 



( 41 ) 

The above differential equation encapsulates the en- 
tropic effect, brownian motions of the monomers and the 
monomer-monomer interaction. We believe that it would 
serve to be a natural platform for the study of single 
polymer dynamics such as in the scenarios of polymer 
collapse JllJl^J^lglLl^j |23|,j24l SH] , translocation 

and adsorption [3(| . 
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